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PRACTICA 1 

Viga de Aluminio en Voladizo 


Description de la practica: 

La practica realizada el lunes 20 de febrero consistio en analizar las caracterlsticas dinamicas de 
una viga de aluminio, se considero la viga como empotrada en un extremo y libre en el otro. 
Para simular esta condicion se utilizo una sargenta para simular el empotramiento en un 
extremo de una mesa y se fijo un sensor de aceleracion en el extremo libre de la viga. El sensor 
estaba conectado por un cable a un acondicionador de la serial de tres canales con filtro paso 
bajo y amplificador, la salida del canal 1 de este conectado a una caja de adquisicion de datos 
de 16 canales que funciona a 16 bits. 

Los objetivos de la practica son: 

Describir los elementos de la instrumentation empleados en la practica e indicar los 
parametros empleados en la realizacion de las medidas. 

Obtener la frecuencia propia de la estructura en el dominio del tiempo y analizar su 
variacion. 

Obtener el amortiguamiento del sistema, por los metodos de Decremento Logaritmico 
y Transformation de Hilbert. 

Representar numericamente el problema para ajustar la propiedad de la rigidez del 
apoyo empleando uno de los metodos explicados en clase (Matricial o Rayleigh) 

1. Description de la Instrumentation. 

Sensor. 

Para la realizacion de la practica se uso como primer elemento un sensor de aceleracion 
o acelerometro de la marca comercial ENDEVCO, modelo 7752-1000. A continuation, se 
presenta una imagen tomada el dia de la practica. 


Imagen 1. Sensor de Aceleracion. 

Modelo: 2252.1000 
S/N: 10401 

Sensibilidad: 1039 mV/g 

Suministro de corriente: 4 mA 

Sensibilidad Transversal: 0.8% 
(max) 
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El sensor usado en la practica es un acelerometro con partes electronicas integradas, 
disenado para bajas frecuencias y bajos niveles de aceleracion, ideal para el uso en 
analisis modal. 

Por su pequeno tamano, este tipo de sensores de aceleracion ubican su rango de 
frecuencias de resonancia en valores muy elevados por lo que incurrir en estos rangos 
en mediciones estructurales referidas a la ingenieria civil resulta muy dificil. 

A continuacion, se presentan las caracteristicas geometricas del sensor: 



Es importante la consideracion de la 
geometria y todas las caracteristicas 
propias de los sensores para no sufrir 
percances relacionados a la 
incompatibilidad de conectores, cables, 
fijacion en la base, etc. 


STANDARD TOLERANCE 
INCHES : MILLIMETERS'' 

XX D3 \X« 

xxx - -, 1 - md : xx - **■ _zs; 


Imagen 2. Caracteristicas geometricas del sensor 
Se resumen las especificaciones tecnicas proporcionadas por el fabricante: 


Caracteristicas Dinamicas. 


Rango de trabajo 

g 

±5 

Sensibilidad ±20% 

mV/g 

1000 

Frecuencia de Resonancia 

kHz 

24 

Sensibilidad Transversal 

% 

8 

Caracteristicas de salida. 

Voltaje de salida (DC output BIAS Voltage) 

Vdc 

+7.5 a + 14 

Impedancia de salida 

0 

<250 

Voltaje de salida 

V 

±5 

Recuperacion de sobrecargas 

psec 

<10 
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Requerimientos de poder. 
Suministro de voltaje 

Vdc 

+18 a +24 

Suministro de corriente 

mA 

+2 a +10 

Tiempo de precalentamiento 

Sec 

< 10 

Proteccion de sobrecarga. 

g 

2000 

(no danaran los componentes internos) 



Caracteristicas ambientales. 



Rango de temperatura 

°C 

-20 a +85 

Vibracion sinusoidal 

g pk 

1000 

Limite de impacto 

g pk 

2000 

Caracteristicas fisicas. 



Peso 

G 

33 

Material de recubrimiento 


Acero inoxidable y titanio 

Conector 


10-32 conector superior anclado en 



carcasa, compatible con cables 



serie 3000 

Torque de montaje 

Nm 
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Acondicionador de la Serial. 

La serial de salida del sensor necesita ser amplificada acondicionada para su 
transmision por cable hasta el punto de recepcion, de esta manera evitar que el 
ruido producido por cualquier elemento externo se adhiera a la senal 
distorsionando la lectura. 

Ademas, es necesaria la introduccion de un filtro de paso bajo o antialiasing, para 
nuestro caso en concreto, se utilizo un acondicionador de la marca EC Systems 
modelo PA3000 que funciona como Acondicionador, Amplificador y Filtro 
simultaneamente. 



Imagen 3. Acondicionador usado en la practica. 

Para nuestra practica utilizamos el filtro de paso bajo para que funcione como un 
filtro antialiasing. 
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Este equipo cuenta con los circuitos de alimentacion necesarios para el correcto 
funcionamiento del sensor. Cada canal independiente permite una seleccion de 
rango de amplificacion de lx, lOx o lOOx y los filtros de 1 khz de paso bajo y 10 Hz 
de paso alto. En la practica trabajamos con amplificacion de 1. 

Especificaciones tecnicas. 

Numero de canales 3 


Caracteristicas de entrada 

IEPE24V DC /2.4 mA 

Conector de entrada / salida 

BNC/BNC 

Impedancia de entrada 

100 kQ 

Impedancia de salida 

1200 0 

Ganancia 

lx, lOx y lOOx 

Error de ganancia 

<0.5% 

Rango de salida 

±10 V 

Distorsion 

<0.1% 

Desfase 

<10 mV 

Rango de frecuencia (- 3 dB) 

0.5 Hz -100 kHz 

Filtro de paso bajo 

1kHz, 12 dB / okt 

Filtro de paso alto 

10 Kz, 12 dB / okt 

Indicadores LED 

On / off 

Abierto / Sensor de circuito corto 
Rango sobrepasado 
Ganancia seleccionada 
Cargando 
Filtros habilitados 

Suministro de poder 

12 V DC/400 mA 


Caja de adquisicion de datos. 

Una caja de adquisicion de datos de la serial de 16 canales, producida por National 
Instruments de 16 bits. Tiene entradas y salidas analogicas. 

Para la realizacion de la practica, se uso el canal 0 para la recogida de datos. 

Cumple la funcion de digitalizar la serial para que pueda ser recibida por el ordenador y 
ser usada para los calculos necesarios a posteriori. 

Para la recogida de datos se uso una aplicacion o software desarrollado por los 
integrantes del equipo tecnico del departamento usando LabVIEW. 

La caja tiene una sensibilidad de 1 mV/g. 
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Imagen 4. Caja de adquisicion de datos. 

2. Frecuencia propia en el dominio del tiempo y variation. 

Datos necesarios para la practical 
Modulo de Young: E=7E 10 N/m 2 
Coeficiente de Poisson: p=0.35 
Densidad: p=2650 kg/m 3 
Datos recogidos en la practica 
Dimensiones de la viga: 

L=0.35 m 
B=0.0512 m 
E=0.0035 m 

Para poder representar la grafica del sistema en vibration libre usamos el software 
MATLAB. Para este fin primero graficamos todos los datos obtenidos, usando el criterio 
de seleccionar la "corneta" mas limpia y que represente mejor el movimiento de la 
estructura depuramos las demas para facilitar la lectura, interpretation y la velocidad 
de procesamiento. 

A continuation, se presenta el codigo utilizado, las graficas representativas y los 
resultados: 

%frecuencia de muestreo 

fs=4096; 

dt=l/fs; 

%en el primer tramo 

nTdl=14.5; %numero de ciclos contados 

Tll=140; %tiempo de inicio de conteo s 

T21=141; %tiempo de fin de conteo s 
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fdl=nTdl/ (T21-T11) ; %frecuencia natural de la estructura Hz 

wdl=2*pi*fdl %frecuencia angular natural en rad/s 

Tl=l/fdl 


%en el segundo tramo 
nTd2=l 4 . 5 ; 

T21=140; 

T22=141; 

f d2=nTd2 / (T22-T21) ; 
wd2=2 *pi* f d2 
T2-l/f d2 


%numero de ciclos contados 
%tiempo de inicio de conteo s 
%tiempo de fin de conteo s 
%frecuencia natural de la estructura Hz 
%frecuencia angular natural en rad/s 
%periodo de la estructura 


%en el tercer tramo 
nTd3=l 4 . 3 ; 
T31=142.5; 
T32=143.5; 
f d3=nTd3/ (T32-T31) ; 
wd3=2 *pi* f d3 
T3— 1/ f d3 


%numero de ciclos contados 
%tiempo de inicio de conteo s 
%tiempo de fin de conteo s 
%frecuencia natural de la estructura Hz 
%frecuencia angular natural en rad/s 
%periodo de la estructura 


%en el cuarto tramo 
nTd4=14 .25; 

T41=144; 

T42=145; 

f d4=nTd4 / (T42-T41) ; 
wd4=2 *pi* f d4 
T4=l/fd4 


%numero de ciclos contados 
%tiempo de inicio de conteo s 
%tiempo de fin de conteo s 
%frecuencia natural de la estructura Hz 
%frecuencia angular natural en rad/s 
%periodo de la estructura 


deltaT= ( (T4-T1) /Tl) *100 %porcentaje de incremento del periodo 



Imagen 5. Grafica del movimiento de la estructura, Aceleracion vs. Tiempo 
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Aceleracion(m/s2) 



Imagen 6. Lineas auxiliares para contar los ciclos en los periodos de tiempo escogidos 



Imagen 7. Lineas auxiliares para contar los ciclos en los periodos de tiempo escogidos 

Los resultados de frecuencias amortiguadas obtenidos fueron: 
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fdl = 14.5 Hz 
wdl = 91.1062 rad/s 
Tl= 0.069 s 

fd2 = 14.5 Hz 
wd2 = 91.1062 rad/s 
T2=0.69 s 

fd3 = 14.3 Hz 
wd3 = 89.8495 rad/s 
T3 = 0.0699 s 

fd4 = 14.25 Hz 
wd4 = 89.5354 rad/s 
T4 = 0.0702 s 

deltaT = 1.7544% 

Podemos observar que a medida que el movimiento de la estructura es amortiguado la 
frecuencia de la misma tiende a bajar. Puede ser una consecuencia de la disipacion de 
la energi'a o consecuencia del rozamiento con el aire. 

De manera inversa, el periodo de la estructura tiene a crecer en una razon de 1.7544%, 
corroborando la no linealidad de los sistemas. 

De la primer grafica tambien podemos observar que se sobrepaso el Ifmite de 
aceleracion seleccionado para el ensayo, esto es evidente al ver la formacion de una 
meseta durante los primeros ciclos. 

3. Amortiguamiento y comparacion con Decremento Logaritmico y transformation de 
Hilbert. 

Por decremento logaritmico: 

Primero se identified los siguientes puntos maximos en la serie, los valores de x 
corresponden a los del tiempo y los de y a la aceleracion. 

En la serie fe figuras presentadas abajo, se muestran los puntos de donde se obtuvieron 
valores de tiempo y aceleracion para aplicar el metodo de decremento logaritmico. 
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3eleracion(nrVs2) Heracion(m/s2) 
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X: 141 .5 
Y: 16.22 


2do punto 






4to punto 


X: 143.5 
Y: 1.334 





Con ayuda de un codigo de Matlab se calculo los valores de amortiguamiento para cada 
uno de los intervalos entre puntos, siendo 5 en total, adicionando un valor de 
amortiguamiento calculado para el intervalo entre T1 y T5. El codigo de Matlab 
representa el calculo por decremento logaritmico. 

x ± = x(ti) = Z) 1 e _ ^ wt 1 

x ± = x(t ± ) = 

y 

In — = -^(ti - t 2 ) ; t 2 = t 1 + nT d 

*2 

Por eso el amortiguamiento: ( = — In — 

2tztl X2 

A continuacion, se presenta el codigo: 
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%Valores maximos de Aceleracion y tiempos correspondientes 


T 1 = 1 4 0 . 4 ; 
a 1=4 4 . 51 ; 
f — 1 4 . 5 ; 


%tiempo del primer punto s 

%aceleracion maxima del segundo punto m/s2 
%frecuencia en Hz 


T2=l 4 1.5; 
a2=l 6.22; 


%tiempo del segundo punto s 

%aceleracion maxima del segundo punto m/s2 


T3=l 42 . 5 ; 
a3=5 .21; 


%tiempo del tercer punto s 

%aceleracion maxima del tercer punto m/s2 


T4=143.5; %tiempo del cuarto punto s 

a4=1.334; %aceleracion maxima del cuarto punto m/s2 


T5 = l 4 7.5; 
a5=0 .1175; 


%tiempo del segundo punto s 

%aceleracion maxima del segundo punto m/s2 


%calculo del amortiguamiento por el decremento logaritmico 

z 1= ( 1 / ( 2 *pi* (T2-T1 ) *f) ) * (log (al/a2 ) ) 

z2= ( 1 / ( 2 *pi* (T3-T2) *f ) ) * (log(a2/a3) ) 

z3= ( 1 / (2*pi* (T4-T3 ) *f ) ) * (log(a3/a4) ) 

z 4= ( 1 / ( 2 *pi* (T5-T4) *f ) ) * (log(a4/a5) ) 

z5= ( 1/ (2*pi* (T5-T1) *f) ) * (log (al/a5 ) ) 

deltaz= ( (z5-zl) /zl) *100 


Donde: (Tj-Ti)*f = al numero de ciclos entre los dos puntos seleccionados. 

Del codigo anterior se obtuvieron los siguientes resultados: 

zl = 0.0101 o 1.01% 
z2 = 0.0125 o 1.25% 
z3 = 0.0150 o 1.50% 
z4 = 0.0067 o 0.67% 
z5 = 0.0092 o 0.92% 
deltaz = -8.8806 

Es evidente que los valores obtenidos de amortiguamiento son muy pequenos, ninguno 
superando el 2% e incluso tendiendo a ser mas pequenos. 

Se puede observar que el amortiguamiento entre intervalos no es el mismo, corroborando la 
no linealidad del sistema, ademas, se puede apreciar que cuanto mas se reduce la cantidad de 
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movimiento del sistema, el amortiguamiento tiende a tomar valores menores como es el valor 
de z4 de 0.67%. 

Cuando se calcula un valor de amortiguamiento usando los valores extremos entre los puntos 
maximos escogidos (de T1 a T5) se obtiene un valor de amortiguamiento menor al de los 
obtenidos cuando la amplitud del movimiento del sistema es mayor. 

El decremento del amortiguamiento entre periodos de tiempo es de -8.8806 %. 

El hecho de que el amortiguamiento del sistema sea tan bajo, corrobora la hipotesis de que la 
frecuencia natural amortiguada es igual a la frecuencia natural de la estructura. 

Por la transformation de Hilbert: 

Para la realizacion de la transformacion de Hilbert usamos un comando de Matlab muy 
sencillo, ademas de eso ajustamos una recta con una correlacion suficiente para conseguir un 
valor de amortiguamiento aceptable. A continuacion, se presenta el codigo empleado: 


%transf ormacion de Hilbert 
%Datos Obtenidos en la Clase 

datos=dlmread ( ’ datos_practica_depurados . txt ' ) ; 

fs=4096; %frecuencia de muestreo 

dt=l/fs; %paso de tiempo 

%datos traspuestos 
t=datos ( : , 1 ) ’ ; 
acel=datos ( : , 2 ) ’ ; 

%transf ormada de Hilbert 
y=hilbert (acel) ; 
yhilb=log (abs (y) ) ; 

tajust=(140.5:dt:147.5) ; 
yajust=yhilb (0.5*fs:7.5*fs) ; 

c=polyf it (tajust^yajust, 1) 

%ajuste de la recta 
tp=linspace (140. 5^147. 5 f 7*fs) ; 
yp=polyval (c, tp) ; 

%grafica de la recta ajustada 

figure (2) 

plot (tajust,yajust) 

hold on 

plot (tp, yp) 

hold on 

c = -0.9631 139.3568 

Usando el valor "a=-9631" correspondiente a la correlacion sustituimos en la formula 
encontrada, calculamos el amortiguamiento. 

zhil=- (-0.9631)/ (2*pi*f ) 
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Obteniendo el siguiente resultado: 

zhil = 0.0106 o zhil = 1.06% 



El resultado del amortiguamiento obtenido por este metodo es muy aproximado a los 
resultados obtenidos por el metodo de decremento logaritmico. 

4. Ajuste de la rigidez usando el metodo de Rayleigh o Matricial 

Para hacer el ajuste de la rigidez se realizo un calculo sencillo por un metodo de prontuario, 
de manera que aproximando la relacion de rigideces del empotramiento con la de la viga se 
puede encontrar el valor de la frecuencia del sistema. 

Entrando con la siguiente relacion en un prontuario de Blevins e interpolando: 

M 

— = 0.1985 

mL 

Suponiendo una relacion de rigidez entre el apoyo y la viga de 10 e interpolando: 

A w = 1.5445 

Suponiendo un empotramiento perfecto: 

hnf = 1-6709 

De la relacion: 


A = 
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Se obtienen como resultados: 


flO = 16.0940 Hz finfi = 18.8360 Hz 

Ambos superiores a la frecuencia obtenida experimentalmente (14.5 Hz), de donde se 
deduce inmediatamente que la relacion entre rigideces es menor a 10. 

Posteriormente se hizo un calculo de la frecuencia por el metodo de Rayleigh, este calculo 
supone condiciones de empotramiento perfecto, por lo que el resultado deberia aproximarse 
bastante al obtenido usando el prontuario donde la relacion de rigideces tiende a infinito. 
Para facilitar el calculo se uso un codigo muy corto de Matlab. 

%calculo de la frecuencia por el metodo de rayleigh 

fray= (1/ (2*pi) ) * ( (pi/ (2*L) ) A 2) * ( (E*I) / ( ( ( (3*pi)-8) /pi) *m+(2/L) *M) ) A 0.5 

que arrojo como resultado un valor de: 

fray = 18.0493 Hz que contrastado con el de finfi = 18.8360 Hz se aproxima bastante bien. 

Luego se uso una programacion proporcionada en Matlab para el calculo de los tres primeros 
modos de vibracion por el metodo matricial. A continuacion, se presenta dicho codigo: 

clear 

clc 


E=7el 0 ; 
b=0 . 0512 ; 
esp=0 . 0035 ; 

L=0 . 35 ; 

1= (b*esp A 3 ) /12 ; 
rho=2 650 ; 
m=b*esp*rho; 

M=0 . 033; 


%N/m2 

%ancho de la viga en m 
%espesor de la viga en m 
%longitud de la viga en m 
%Inercia 

%densidad del aluminio, kg/m3 
%masa distribuida de la viga kg/m 
%masa del sensor kg 
klO = 10*E*I/L; % N-m(rigidez del apoyo) 
kinfi = 10 A 10*E*I/L; % N-m(rigidez del apoyo) 
ndiv= 20; % divisiones de la viga 
Le = L/ndiv; % matrices de rigidez elementales 
Ke = E*I/Le A 3* [ 12 6*Le -12 6*Le 
6*Le 4*Le A 2 -6*Le 2*Le A 2 
-12 - 6*Le 12 - 6*Le 
6*Le 2*Le A 2 -6*Le 4*Le A 2]; 

Me = m^Le/ 420* [ 156 22*Le 54 -13^Le 
22 *Le 4^Le A 2 13^Le -3^Le A 2 
54 13^Le 156 -22^Le 
-13^Le -3^Le A 2 -22*Le 4*Le A 2]; 


KK = zeros ( 2 * (ndiv+1 ) , 2 * (ndiv+1 ) ) ; MM = zeros ( 2 * (ndiv+1 ) r 2 * (ndiv+1 ) ) ; 
KK(2 a 2) = klO; % rigidez concentrada apoyo 
MM (end-1 , end-1 ) = M; % masa concentrada 

for i=l:l:ndiv % ensambla matriz de masa y rigidez 
idx = ( ( i — 1 ) * 2 + 1 : ( i — 1 ) *2 + 4) ; 

KK(idx,idx) = KK (idx, idx) +Ke ( 1 : 4 , 1 : 4 ) ; 

MM (idx, idx) = MM ( idx, idx) +Me ( 1 : 4 , 1 : 4 ) ; 
end 

T=eye (2* (ndiv+1 ), 2* (ndiv+1 )); T (:, 1 )=[] ; % gdlactivos 
K=T'*KK*T; M=T'^MM^T; % condiciones de contorno 
[V,D] = eig(K,M); % calculo de autovalores 
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frq= diag(D) . A .5/2/pi;% frecuencias en Hz 
mod=T*V; % expansion gdlen modos 

% representa 3 primeros modos 

figure(l); plot ( ( 0 : 1 : ndiv) ' *Le , mod ( ( 1 : 2 : 2 * (ndiv+1 ) ) , 1 : 3 ) ) ; 
legend (num2str ( f rq ( 1 : 3 ) ) ) ; 

Obteniendo una grafica que representa los tres primeros modos de vibracion de la estructura: 
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Imagen 8. Tres primeros modos de vibracion. 

De la grafica se puede constatar que la frecuencia correspondiente al primer modo de vibracion 
esta aproximada a la frecuencia natural obtenida experimentalmente, siendo un 1.045% mayor. 

Finalmente, para determinar la rigidez del apoyo, se volvio a la relacion proporcionada por 
Bleming: 


A = 


% 

2 * nL 2 



Donde conocido el valor de / x (obtenido experimentalmente) se despeja la variable A. 


Obteniendo un valor de A = 1.466. De manera que interpolando en filas y columnas se obtiene 
una relacion entre rigideces de: 


kL 

- = 7.5456 


Despejando el valor de rigidez del apoyo se obtiene un valor de: 


k = 276.0684 N/m 

Resultado en un valor considerablemente bajo para ser supuesto como un empotramiento 
perfecto. 
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